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ABSTRACT 

We estimate the Rees-Sciama (RS) effect of super structures on the cosmic microwave background 
(CMB) temperature fluctuations and identify a related effect on galaxy redshifts. By numerically 
solving the geodesic equation, we find that both superclusters and supervoids can decrease the tem- 
perature of the CMB by several micro Kelvin in the central region and increase the temperature 
slightly in the surrounding area due to the RS effect. The two components of the RS effect, redshift 
and gravitational time delay, largely cancel each other, leaving an equivalent but much smaller effect 
on the CMB photons that started out at the same time from the distorted last scattering surface. For 
galaxies, the time delay effect is separable from the redshift effect, and the slight change to the redshift 
induced by super structures can be at the percent level of large-scale rms bulk velocities, which might 
only be detected statistically. On much smaller scales, a tiny redshift difference between two images 
of a strongly lensed source should exist in general, which is related to the Hubble expansion rate at 
the source redshift. However, as Loeb (1998) pointed out, observational issues and the proper motion 
of the structure would make such a measurement impossible. 

Subject headings: cosmic microwave background — cosmology: theory - galaxies: distance and red- 
shifts — large-scale structure of universe — relativity 



1. INTRODUCTION 

Measurements of the cosmic microwave background 
(CMB) with the Wilkinson Microwave Anisotropy 
Probe ( WMAP, Bennett et al. 2003) have played 
a crucial role in establishing the concordance cos- 
mological model, ACDM - a flat eold-dark-matter 
universe with roughly three quarters of its content 
in the form of dark energy (e.g., Spergel et al. 2003). 
Besides the primary anisotropics, one also expects 
secondary CMB anisotropics arising from late-time 
gravitational effects, such as the integrated Sachs-Wolfe 
(ISW) effect (Sachs & Wolfe 1967), and scattering 
processes, such as the thermal Sunyacv-Zcl'dovich 
effect (Sunyaev & Zeldovich 1972). Detections of 
these effects (e.g., Birkinshaw, Hughes, & Arnaud 
1991; Carlstrom, Joy, & Grego 1996; 

Fosalba, Gaztahaga, & Castander 2003; 
Boughn & Crittenden 2004; Afshordi, Loh, & Strauss 
2004; Padmanabhan et al. 2005; Giannantonio et al. 
2008; Granett, Neyrinck, & Szapudi 2009) are useful 
confirmations of the ACDM model. 

The ISW effect on CMB photons is caused by the time 
variation of gravitational potentials as the photons travel 
through them. In order to have a measurable ISW effect, 
the photon travel time must be an appreciable fraction 
of the time scale, over which the potentials vary signif- 
icantly. On large scales, overdensitics in an Einstein-de 
Sitter (EdS) universe grows at the same rate as the cos- 
mic expansion, so that the linear ISW effect vanishes. 
However, in a ACDM universe, the cosmic expansion rate 
exceeds the linear growth rate, causing a decay of poten- 
tials. A photon traveling through a decaying potential 
well (wall) gains (loses) energy. Therefore, detection of 
the linear ISW effect provides a piece of evidence for dark 
energy in a flat universe. 



Nonlinear clustering causes extra evolution of the po- 
tentials against the background expansion, so that the 
ISW effect is present even in an EdS universe. This 
effect is also known as the Rees-Sciama (RS) effect 
(Rees & Sciama 1968). On large scales, the nonlinear 
contribution (or the intrinsic RS effect) to the full ISW 
effect is expected to be subdominant to the linear ISW 
effect in a cosmological constant (or dark energy) domi- 
nated universe (Seljak 1996; Tuluie, Laguna, & Anninos 
1996). 

In a recent paper, Granett, Neyrinck, & Szapudi 
(2008) report a mean CMB temperature change of 
| AT| = 9.6 ± 2.2 /iK caused by superclusters and super- 
voids at z ~ 0.5 (hereafter, we drop the prefix "super" 
where there is no ambiguity). These structures have radii 
of ~ 100/i _1 Mpc or 4° at z ~ 0.5. The scales involved 
are large enough, so that the linear ISW effect should 
be the major source of the temperature change in the 
ACDM model. However, the result is somewhat puz- 
zling, because the temperature change is twice of what is 
estimated from simulations (Granett et al. 2008). Fur- 
thermore, there is no significant difference between the 
mean temperature of the clusters and that of voids in the 
reconstructed ISW map (Granett, Neyrinck, & Szapudi 
2009). 

Given the studies above, it is of interest to reexamine 
the different mechanisms, through which cosmic struc- 
tures alter the CMB temperature. We focus on the 
intrinsic RS effect in this paper. Our goal is not to 
reconcile the reported mean CMB temperature change 
with theory, but we note that the thermal Sunyaev- 
Zel'dovich effect and the proper motion effect of struc- 
tures (Birkinshaw & Gull 1983; Gurvits & Mitrofanov 
1986; Stebbins 1988) are not likely to explain the discrep- 
ancy. In WMAP bands, unresolved thermal Sunyaev- 
Zel'dovich signals would reduce the CMB temperature 
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of clusters with respect to that of voids and hence sup- 
press \AT\/T, whereas the proper motion of a cluster or 
void would produce a dipole pattern without affecting 
the mean temperature. 

We adopt the EdS universe as the background cosmol- 
ogy to isolate the intrinsic RS effect (we drop the word 
"intrinsic" hereafter). Because the growth rate of struc- 
tures varies with the assumed cosmology, the RS effect 
is model dependent. For example, the combined RS and 
proper motion effects decrease by a factor of ~ 2 from the 
EdS universe to an open universe with matter fraction 
0.3 (Tuluie et al. 1996). 

One approach to investigate the RS effect is to analyze 
its power spectrum (e.g., Seljak 1996). The linear ISW 
effect increases the CMB temperature fluctuation and 
hence the power of the modes by raising the temperature 
of overdense regions and lowering that of underdense re- 
gions. Because the RS effect produces a net temperature 
decrement for both clusters and voids (see Section 5 for 
an explanation), it would not increase the temperature 
fluctuation as much even if it had the same amplitude 
as the linear ISW effect. Thus, power spectrum analyses 
may underestimate the importance of the RS effect. 

Ray-tracing through iV-body simulations can pro- 
vide a realistic estimate of the full ISW effect (e.g., 
Tuluie et al. 1996; Maturi et al. 2007; Cai et al. 2009, 
2010; Smith et al. 2009) To do so exactly, one needs to 
solve for the metric and its spatial-temporal derivatives, 
which appear in the geodesic equation. The ISW effect 
is roughly \ AT\/T = \ Az\/(l + z) ~ lO^-lO" 5 for super 
structures, so that a redshift error of the order 10 ~ 7 or 
larger accumulated along the photon geodesic, or a com- 
parable error from approximating the photon geodesic 
with a straight coordinate ray, can be a significant con- 
tamination to the results. Therefore, it is computation- 
ally challenging to trace photon geodesies in TV-body sim- 
ulations. 

In this paper, we study the effect of individual struc- 
tures with toy models (see, e.g., Thompson & Vishniac 
1987; Martinez-Gonzalez & Sanz 1990). The advantage 
of this approach is that (1) for certain class of mod- 
els, one can obtain analytic solutions of the underly- 
ing metric without worrying about spatial, temporal, 
or mass resolutions and (2) in models with symme- 
try, certain null geodesies can be computed in differ- 
ent ways, so that one can examine the precision of 
the calculations. Specifically, we model the clusters 
and voids with the spherically symmetric Lemaitre- 
Tolman-Bondi (LTB, Lemaitre 1933; Tolman 1934; 
Bondi 1947, referenced herein) solution and solve generic 
photon geodesies through these structures (see also 
Panek 1992; Alnes & Amarzguioui 2006; Marra et al. 
2007; Valkenburg 2009). 

Separately, the RS effect has been studied for 
compensated voids and clusters (Inoue & Silk 2006; 
Tomita & Inoue 2008; Sakai & Inoue 2008) and for com- 
pensated shells (Afshordi, Slosar, & Wang 2011) in the 
ACDM universe with thin-shell approximation, pertur- 
bative calculations, and numerical calculations. The fo- 
cus of this paper is on the effect caused solely by the evo- 
lution of the structures, so we do not include the cosmo- 
logical constant A except for a simple case in Section 6. 
Our results are in qualitative agreement with those in 



the aforementioned works, and we extend the study to 
uncompensated structures as well. 

As Rees & Sciama (1968) pointed out, redshift and 
time delay are the two major components of the RS ef- 
fect. The redshift component includes both compensa- 
tion for the time delay component and the evolution of 
the potential. Since the effect of the potential evolution 
is usually subdominant, the two components of the RS 
effect often shift the CMB temperature with comparable 
magnitudes but in opposite directions. 

Galaxy redshifts are directly measurable, so fractional 
changes to galaxy redshifts due to intervening structures 
are larger than that to the CMB temperature. This 
foreground-induced change of galaxy redshifts is closely 
related to but different from the RS effect or the ISW 
effect in general because of the separation of the red- 
shift effect and the time delay effect on galaxies. How- 
ever, unlike the CMB, which has a standard temperature, 
galaxies are spread over redshift space. Detecting such a 
redshift change is far more difficult than that of the RS 
effect. 

Images of a strongly lensed source should have slightly 
different redshifts in general. Despite the extreme diffi- 
culty in measuring such redshift differences (Loeb 1998), 
we find in the ideal case that the redshift difference be- 
tween two images of the source is considerably larger 
than that between two epochs of observations of the same 
image separated by the amount of the time delay between 
the two images. The latter is known as the Sandage-Loeb 
test (Sandage 1962; Loeb 1998). 

The rest of the paper is organized as follows. Section 2 
provides a brief introduction to the LTB solution and de- 
scribes models of clusters and voids. Section 3 gives the 
details of solving the null geodesies numerically. Results 
of the RS effect on CMB temperature profiles and tiny 
perturbations to galaxy redshifts caused by intervening 
structures are presented in Sections 4 and 5, respectively. 
Further discussion is made in Section 6. 

2. LTB MODELS OF SUPER STRUCTURES 

We model the structures as spherically symmetric, 
dust-filled objects embedded in the EdS universe. The 
line element in such models is described by the LTB met- 
ric 

R ,2 (t,r)dr 2 



ds' 



-dt l 



1 - K(r)r 2 



R 2 (t,r)dtt 2 



(1) 



where R(t, r) is the angular diameter distance of the co- 
ordinate r as viewed from the center, K(r) is the curva- 
ture function, a prime denotes a partial derivative with 
respect to r, and the speed of light has been set to unity. 
The evolution of R(t,r) is determined by K(r) and the 
mass function M(r) 



R 2 (t,r) 



2GM(r) 
R(t,r) 



K{r)i 



(2) 



where G is Newton's constant, and an overdot stands for 
a partial derivative with respect to t (hereafter, we sup- 
press the variables t and r if there is no ambiguity in the 
context). The mass function is related to the acceleration 
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and the coordinate density is given by 

w 



4TrR 2 R' 



(4) 



The general parametric solution of Equation (2) for 
K > is 



R 



t-t B (r) 



GM 
GM 



{u — sin it). 



(5) 
(6) 



The bigbang time, ig(r), defines a coordinate surface, 
R[tB(r), r] = 0, at bigbang and is determined by M and 
K up to a scaling of r. It ensures that every part of the 
inhomogeneous universe contracts to a singular point at 
the same coordinate time as we look back. The solution 
for K < is obtained by replacing u with iu. In the 
special case, K = 0, we have 



9GM \ 



1/3 



[t-t B {r)f/\ 



(7) 



If both K and M are constant, i.e., the universe is homo- 
geneous and isotropic, the solution has the form R = ar 
with a being the scale factor of the cosmic expansion. 
More specifically, we reserve a for the scale factor of the 
EdS universe, i.e., the background, in this paper. 

A particular realization of the LTB model is gener- 
ated by specifying K (or M) and initial values of R and 

R. Following Marra et al. (2007) and Paranjape & Singh 
(2008), we set 



R{U,r) = a*r and R(U,r) = aiHir, 



(8) 



where a» = a(ij) and Hi = H(ti) are, respectively, 
the scale factor and Hubble expansion rate of the EdS 
background at the initial time ti. For convenience, we 

choose 014 = 10~ 3 and £Z,a^ 2 = Hq = 50kms _1 Mpc -1 . 
This rather small value of the Hubble constant is of- 
ten adopted in the EdS universe to satisfy the age con- 
straint, and, for this work, it also roughly matches the 
angular size of the model structures in the EdS uni- 
verse with that in Granett et al. 2008 for the same linear 
size. Because we set the linear size of the structures 
in units of Mpc, a lower H would make these struc- 
tures smaller relative to the observable universe and re- 
duce their RS signal. Results with a larger Hq are given 
in Section 4 for comparison. Equation (8) ensures that 
the overdense (underdense) region always contracts (ex- 
pands) against the background after ti (but before the big 
crunch, if applicable) and thus becomes more and more 
overdense (underdense). One may also set up voids with 
the overdensity in underdense region not always decreas- 
ing with time (e.g., Alnes, Amarzguioui, & Gr0n 2006; 
Kolb, Marra, & Matarrese 2008). 

We study both compensated and uncompensated LTB 
models. The former requires the mean density of the 
region to match that of the EdS background beyond a 
finite radius, and hence there is no effect on the rest 
of the universe due to Birkhoff's theorem. We design 
our toy models to have a smooth density profile (i.e., at 
least continuous in its first derivative with respect to r) 
with a uniform inner region (region I) inside a radius n, 
so that numerical results there can be checked against 



analytic solutions. The compensated model is realized 
by setting K = beyond an outer radius, ni- We refer 
to the region between r\ and rn as region II and that 
outside rn as region III, which evolves exactly as the 
EdS universe. For the uncompensated model, we set the 
initial coordinate density in region III to the initial EdS 
background density pEds,i- 

Because the geodesic equation involves R", a discontin- 
uous R" (such as those conforming with top-hat density 
profiles) will deflect photons in coordinate space (but not 
in physical space) as they travel across the discontinuity. 
This does not have theoretical impact on the investiga- 
tion, but to ensure the numerical precision, one must ac- 
curately determine the 4-coordinate x^ where the photon 
crosses the discontinuity and then adjust its wavevector 
fc M , so that its physical velocity remains the same and 
fc M fc M = across the discontinuity. We avoid such a nu- 
merical issue by designing a curvature function that is 
continuous to at least second order (so are R and M), 
which renders smooth density profiles via Equation (4). 
The functional forms of K and corresponding density 
profiles are given in Appendix A. In summary, the pa- 
rameters of the toy models are r\, rn, and the initial 
overdensity of region I, 5 Ij4 = pu//9Eds,; - 1- 

Although there is no peculiarity analytically as the cur- 
vature function K approaches 0, e.g., just inside of ni in 
the compensated models, Equations (5) and (6) are not 
suitable for numerical evaluation when \K\ is very small. 
Observing that q = u|isT| -1 / 2 remains finite as \K\ — > 0, 
we Taylor expand Equations (5) and (6) around u = 
and then replace u with g |7iT | x / 2 when K ~ 0. 

3. GEODESICS 
The geodesic equation is given by 

d 2 x^ „„ dx a dx 



dv 2 



dv dv 



0. 



(9) 



/here 



(t,r,0,(j>), v is the affine parameter, and 



is the Christoffel symbol. In the LTB metric, Equa- 
tion (9) becomes 



dk l 
dv 
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dv 
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dv 
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dv 
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K'r 2 



2(1 - Kr 2 ) 



R' 



2 dR 
~R~db 
2 dR 
~R~db 



k e + sin cos (k?)' 
k^ - 2cot0fcV, 



(10) 

(11) 
(12) 
(13) 



where k^ = dx^/dv is the wavevector, (k n ) = (k e ) + 

sin 2 6 (fc*) 2 , and dR/dv = i?fc t + J R'fc r . We enumerate the 
index \i explicitly as t, r, 8, and 4> instead of numbers for 
clarity. Because of spherical symmetry, we only need to 
study geodesies in the equatorial plane. 

The frequency of a photon with wavevector k 11 as mea- 
sured by an observer (or source) with 4-velocity is 



v cx -k»U u 



(14) 
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If both the observer and the source are comoving in the 
LTB metric, then the photon redshift is given by 



1 



k (t e , r e ) 
fe* {t , r Q ) 



(15) 



where the subscripts e and o denote the event of emis- 
sion and that of observation, respectively, and fc* = 
/c*(i e ,r c ) = 1 is set arbitrarily as the initial condi- 
tion. In actual calculations, we propagate photons back- 
ward from the observer to the source surface by setting 
v —> —v (hence fc M — s- — fc M ) and the initial condition 
fc* = fc*(£ ,r ) = —1. While Equations (9-15) remain 
the same under the reversal of u, others may change sign. 
To avoid confusion, wc only refer to the forward case in 
all the equations and discussions below. 

The RS effect under investigation is of the order 
\AT\/T = \ Az\/(l + z) ~ lCT 6 -lCr 5 , so one must ensure 
that numerical errors in k l is much less than one part in 
a million. Numerical results can be easily checked in 
uniform regions where alternative solutions exist. In ad- 
dition, we perform several general tests that may detect 
numerical errors, which are described as follows. 

By definition, null geodesies obey k^k^ = 0. This is a 
redundant constraint once the initial condition is set, as 
in principle the geodesic equation does not induce viola- 
tion of the condition. However, numerical errors could 
be accumulated. While enforcing redundant constraints 
numerically is a subject of research itself, we simply ad- 
just the time steps so that \k^k^ /ktk 1 ] < 10 -8 when the 
photons exit the systems. This condition is not suffi- 
cient to validate the results, but violation of it indicates 
significant numerical errors. 

In spherically symmetric systems, each geodesic re- 
mains in a plane. One can always rotate the coordinates 
to place the geodesic in the equatorial plane, in which 
case R 2 k^ is conserved. Similarly, R 2 k is conserved if 
k^ = 0. We do not apply the known solution to reduce 
the dimensions of the system. Rather, we use them to 
check the precision of the numerical solutions. The max- 
imum fractional error of the conserved quantity along 
off-center geodesies is found to be ~ 10 -9 . 

Finally, radial geodesies can be calculated directly from 
the LTB metric Equation (1) 



dt 
dr 



= ± 



R' 2 



1 - Kr 2 ' 



(16) 



where outward (inward) geodesies take the positive (neg- 
ative) sign. By differentiating Equation (15) with respect 
to r(v), one gets the evolution of redshift 



dln(l + z) 
dr 



1 dkl 
/c* dr 
1 



1 dv 
fc* dr 



-T* k a kP 1- 

KK L a pk k -± v - 



K 
dv 

R' 



Kr 2 



(17) 



where the sign on the far right side is the same as that in 
Equation (16) (see Mustapha et al. 1998; Celerier 2000, 
for an alternative derivation). Fractional differences be- 
tween numerical solutions of z(r) and t(r) using Equa- 
tions (10-11) and those using Equations (16) and (17) 
are < 10~ 8 when photons exit the model structures. 



4. REES-SCIAMA EFFECT OF SUPER STRUCTURES 

The clusters studied in Granett et al. 2008 have radii 
of roughly 100 /i _1 Mpc and a mean galaxy overdensity of 
0.72 within the overdense region, corresponding to a mat- 
ter overdensity of 0.36 with a galaxy clustering bias of 
roughly 2 for luminous red galaxies (Padmanabhan et al. 
2007; Blake et al. 2008). Thus, our simple LTB models 
of these clusters have an inner radius r% — 140 Mpc 1 , an 
outer radius r\\ = 280 Mpc (205 Mpc) for the compen- 
sated (uncompensated) profile, and a uniform overden- 
sity Si = 0.36 in region I at z = 0.5. For voids, the only 
difference is 5\ = —0.19. Appendix A provides more de- 
tails about these profiles and their corresponding metric 
quantity R. 

To calculate the RS effect, we place a model structure 
between an observer and a source surface and propagate 
photons from the latter to the former (as mentioned in 
Section 3, they are actually propagated backwards). The 
source surface is spherical around the observer in coordi- 
nate space, so that in absence of the structure, photons 
originated from the surface at the same time would reach 
the observer at the same time with the same redshift. 
With the intervening structure though, the photons re- 
ceived at the same time would have started from the 
source surface at different times because of gravitational 
time delay, and their redshifts would have to compen- 
sate for the time delay plus additional changes because 
of nonlinear evolution of the structure. The net change 
in temperature is a sum of the two components: 



AT = —Tomb 



H s At K 



Az 



l + z s 



(18) 



where Tomb = 2.73 K is the present CMB temperature, 
At s is the time difference at the source, H s = H(z s ) is 
the Hubble parameter at the source redshift z s , and Az s 
is the redshift difference "measured" by an observer at 
present. Note that the two components in Equation (18) 
cannot be measured separately in the case of the CMB. 
The time component accounts for the decrease of the 
CMB temperature with time, and Az s bears the sub- 
script s because it varies with the source location. 

In the uniform background universe, a small time dif- 
ference of At translates to a redshift difference of 

Az = -(1 + z)H(z)At. 

Hence, in region III of compensated models, one can sep- 
arate z s into a time-compensation term 

z\ = -(l + z s )H s At s 

and an evolution term z| = Az s — z*. Since z l s exactly 
cancels H s At s in region III, Equation (18) becomes 

Azt 



AT = —Ti 



CMB 



1 



(19) 



In the compensated models, the only place where CMB 
photons are affected by the structure is within r\\ . Hence, 
the RS effect (and ISW effect in general) does not depend 
on the location of the hypothetical source surface as long 
as the compensated structure is enclosed between the ob- 
server and the source. In other words, z| scales linearly 

We adopt h = 0.71 here only to obtain the size of the observed 
structures in units of Mpc. 
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Fig. 1. — CMB temperature shifts due to intervening model struc- 
tures at 2 = 0.5, which have similar sizes (n = 140 Mpc or ~ 4°) 
and overdensitics (<5j = 0.36, —0.19) as the clusters and voids in 
Granett et al. (2008). The average densities of the compensated 
cluster (solid line) and void (dotted line) within the outer radius 
rji = 280 Mpc equal the background density. The uncompensated 
cluster (dot-dashed line) and void (dashed line) have very simi- 
lar inner profiles (r < 1.4rj) as their compensated counterparts 
(see Figure 4). The most prominent feature of the uncompensated 
models is a CMB dipole, which has been removed. 

with 1 + z s in the uniform background. For the uncom- 
pensated models, the density evolution just beyond ru is 
still significantly different from that in the EdS universe, 
so we place the source far away from the model structure. 

Figure 1 illustrates the RS effect of four LTB model 
structures. The main feature is a temperature decrement 
of several fiK in the inner region and a slight temperature 
increase at the outskirt, regardless whether the structure 
is a cluster or a void. This is distinct from the linear 
ISW effect in a universe dominated by dark energy, which 
generally causes a temperature increase for clusters and 
a decrement for voids. The RS effect of the voids may 
seem counterintuitive at first, and we find that it is the 
result of undcrdensitics generally evolving slower than 
the expansion rate in the EdS universe. 

Because a distant comoving observer or emitter in the 
uncompensated LTB models would be seen to have a pe- 
culiar velocity in the EdS background, the most promi- 
nent feature of these models is a CMB temperature 
dipole. Once the dipole is subtracted, the temperature 
profiles are similar to those of the compensated models. 

As mentioned in Section 2, we set a low value for the 
Hubble constant in the EdS background. If we used in- 
stead Ho = 71 kms -1 Mpc -1 , the central temperature 
decrement of the compensated cluster (void) in Figure 1 
would become 9.9/iK (4.2/iK), nearly tripling the decre- 
ment of 3.6/iK (1.5/zK) with H = 50kms _1 Mpc -1 . 

Figure 2 shows the dependence of the RS effect on the 
model parameters n, rn, and S\. We use the difference 
between the characteristic peak and trough of the tem- 
perature profile (T + — T_) in Figure 1 to represent the 
amplitude of the RS effect and show results for compen- 
sated clusters. Because it takes more time for CMB pho- 
tons to go through a larger structure, its potential can 
evolve more to produce stronger RS effects. Structures 
with larger |<5i| evolve faster and also produce larger RS 
effects. These are indeed seen in Figure 2. Moreover, the 
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Fig. 2. — Amplitude of the RS effect of compensated clusters 
as a function of the inner radius with rjj = 2r\ (solid lines) and 
rn = 3n (dashed lines). Thick lines are for clusters with Si = 0.5 
at z = 0.5, while thin lines are for <5i = 0.2. 

RS effect increases fairly rapidly with n, ru, and 6\. 

The several fiK level RS effect in Figure 1 is significant 
compared to the detected CMB temperature shift of ~ 
±10 fiK due to super structures and is comparable to the 
expected maximum linear ISW effect of 4.2 fiK within 
an aperture of 100 /i -1 Mpc in the Millennium simulation 
(Granett et al. 2008). Although the RS effect could bring 
the expected total ISW effect of the supervoids more in 
line with the observed value, it would also widen the 
discrepancy between the observation and the expected 
total ISW effect for superclusters. 

5. EFFECT OF EVOLVING STRUCTURES ON GALAXY 
REDSHIFTS 

Figure 3 illustrates the time delay component (dotted 
lines) and the redshift component (dashed line) of the 
RS effect (circles) as given in Equation (18). These two 
effects are opposite to each other and are of the same 
order, so that the net result is much smaller than cither 
one. In other words, the gravitational time delay in the 
potential is a more dominant effect than the potential 
evolution for super structures. 

Treating recombination as an instantaneous event, 
which suffices our purpose, one can view the RS effect 
and the ISW effect in general as a single redshift effect 
without time delay. This view is fully equivalent to that 
of Rees & Sciama (1968). However, one is no longer com- 
paring photons that have started from the same spatial 
coordinates at different times; rather, the CMB photons 
received at the same time were emitted at the same in- 
stance of recombination, but they must have started from 
different locations. In other words, the last scattering 
surface as viewed by the observer is no longer spherical 
in presence of inhomogeneities. 

To cross-check the results, we trace photon geodesies 
to the same time surface rather than the same spatial co- 
ordinate surface described in Section 4. This is analogous 
to theoretical calculations of the ISW effect by integrat- 
ing over (conformal) time. The resulting change in the 
temperature is then given solely by the redshift differ- 
ence, i.e., a change in k 1 . As shown in Figure 3, this 
single redshift effect (solid lines) gives exactly the same 
RS effect as that from the combination of the time delay 



6 



Zhan 



3. 



60 
40 
20 


-20 
-40 
-60 



_ I I I 1 I I I 1 I I I 1 I _ 

" — .. Compensated cluster ~ 
AT t =-2AtT CMB /3t ~- 


_ , , , | , , , | , , , | , 

Compensated void _ 

~~ ~~ - ^ AT 
_ ^ z 

- \ 

— 

— , 


~~ / ~ 
/ 

- / ' ^T z = — Az T CMB / ( 1 + z) J 
' , , , 1 , , , 1 , , , 1 , " 


AT. 

— o 5x(AT t +AT z ) (same coor.) 

5xAT (same t) ~~ 

, , , 1 , , , 1 , , , 1 , 



> 



a 



2 4 
8 (degree) 



2 4 6 

6 (degree) 



We can estimate the very effect along radial geodesies 
using Equations (16) and (17). The required R and R! 
evolve according to 



Fig. 3. — Two equivalent views of the RS effect: a combination (circles) of the redshift (dashed lines, AT Z ) and time delay (dotted line, 
ATi) effects on photons started at the same coordinate distance from the observer and a single redshift effect of photons started out at 
the same time (solid lines), e.g., the instance of recombination. The left panel and the right panel are for the compensated cluster and the 
compensated void, respectively. The combined RS results (circles) are the same as those in Figure 1 but amplified by 5 times for better 
viewing. For galaxies, their redshift changes are given by (1 + z)AT z /TcmBi i- e -> slightly more redshifted behind clusters and less behind 
voids in the EdS universe. 

and redshift effects (circles). 

Unlike CMB photons, light from a comoving galaxy 
must be emitted from the same coordinates whether 
there is an intervening structure or not. Furthermore, 
galaxy redshifts can be directly measured through spec- 
troscopy, so that the redshift effect is no longer mixed 
with the time delay effect. Without the large cancel- 
lation from time delay, changes of galaxy redshifts can 
be roughly an order of magnitude greater than frac- 
tional changes in the CMB temperature. For the clusters 
(voids) studied in Section 4, the redshift (blueshift) effect 
is merely several kilometers per second, or at the percent 
level of the large-scale rms bulk velocity, estimated to be 
a few hundred kilometers per second (e.g., Zhan & Fang 
2002). Nevertheless, it might be detectable statistically 
in the future and could be somewhat significant to large- 
scale redshift surveys that aim to precisely measure the 
redshift distortion effect. 



A = H$(l-tl m )A- 

AA' K' /Oj 
~A~ + 2i 



2A 2 



A' 



(20) 
(21) 



6. DISCUSSION AND CONCLUSIONS 

We have studied the intrinsic RS effect of super struc- 
tures due to their weakly nonlinear evolution using ide- 
alized LTB models and find that clusters and voids re- 
sembling those in Granett et al. (2008) can cause CMB 
temperature decrements of a few [iK in the EdS uni- 
verse. The amplitude of the RS effect is fairly significant 
compared to the measured CMB temperature changes, 
so one needs to take into account both the RS effect and 
the linear ISW effect when interpreting the observations. 

It is interesting to note that the RS effect is a cen- 
tral CMB temperature decrement, regardless whether 
the structure is a cluster or a void. The reason is that 
generally overdensities grow faster than the cosmic ex- 
pansion in the EdS universe, while underdensities evolves 
slower than the cosmic expansion. 

In the concordant ACDM universe, the linear ISW ef- 
fect is expected to be stronger than the RS effect at low 
redshift when dark energy is dominant. A full calcula- 
tion of the total ISW effect with Equation (9) is rather 
involved, as the metric function R and its first and second 
derivatives are no longer analytic in the ACDM universe. 



where A = r R, and fl m = 0.3 is the present matter 
fraction of the background flat ACDM universe. The ini- 
tial conditions arc given by Equation (8), and the Hub- 
ble constant is set to 70kms _1 Mpc _1 in this case. For 
the compensated cluster (void) in Figure 3, we find that 
AT Z = -21.5 /iK (13.3 /LtK), AT t = 26.6 /iK (-18.7 /M), 
and the net temperature shift AT = 5.1 /iK (— 5.4/iK) 
through its center. These results agree with the gen- 
eral expectation that photons gain (lose) energy going 
through a decaying potential well (wall) in the acceler- 
ating ACDM universe. 

For galaxies, the redshift effect is separable from the 
time delay effect, and the magnitude of the change in 
redshift is much larger than the fractional change in the 
CMB temperature. This redshift effect is due in most 
part to the compensation for the time delay effect, e.g., 
a galaxy seen at an earlier time also has a higher redshift. 
Evolution of the potential plays a minor role in the case of 
super structures. Hence, the effect on the galaxy redshift 
is different from the RS effect or the ISW effect. 

The compensated cluster (void) in Figure 3 changes 
the redshift of a z s ~ 0.6 background galaxy by Az s = 
2.8 x 10~ 5 (-1.7 x 10~ 5 ) in the EdS universe and by 
Az s = 1.2 x 10~ 5 (-7.5 x 10~ 6 ) in the ACDM universe. 
Although measuring such redshift differences is not com- 
pletely out of technical reach, the fact that there is not a 
standard redshift to compare with makes it impossible to 
measure this effect on individual galaxies. Since it is a 
systematic perturbation to galaxy redshifts, one might 
be able to detect the effect by comparing the galaxy 
redshift distribution with and without intervening struc- 
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tures. However, many issues come into play. For ex- 
ample, magnification by the very foreground structures 
also affects how background galaxies are selected into the 
sample and thus change the galaxy rcdshift distribution 
perhaps more significantly. 

While the absolute time delay due to foreground struc- 
tures is not directly observable, one can obtain the rel- 
ative time delay between images of a strongly lensed 
source either by actual measurements with a variable 
source or by estimation with a lens model. This could 
have allowed a short cut to the Sandage-Loeb test that 
has to run for decades. Assuming for the moment that 
the potential does not evolve, i.e., AT = in Equa- 
tion (18), we have 



Av 



cAz^ . cAto 

= -cH a M s = - z -H s , 



1 



(22) 



where Atg is the relative time delay observed at z = 0. 
For a source at z s — 1, the velocity difference correspond- 
ing to a time delay of one year is — 1.8cms _1 . The evo- 
lution of the lens potential also contributes to the red- 
shift difference. Because lenses are more or less virialized 
structures, one might approximate them as static objects 
in physical space and determine their potential evolution 
in the expanding universe. One might also stack CMB 
measurements behind a large number of similar lenses to 
determine the potential evolution. 

For comparison, the Sandage-Loeb test gives (Loeb 
1998) 

1 [H a -H {l + z)}. (23) 



Av 



l + z s 



Here Atg is the time between two epochs of measure- 
ments. In the ACDM universe, H s < Hq(1 + z) at 
z s < 2.5, so the direction of the redshift change in 
the strong-lcnsing case is opposite to the Sandage-Loeb 
test at z s < 2.5. The amplitude in the strong- lensing 
case is almost an order of magnitude larger than that 
in the Sandage-Loeb test, making the former compara- 
tively easier to measure for the same Ato. Moreover, the 
Sandage-Loeb test requires extremely high-precision ab- 
solute calibration of the instruments over a few decades, 
which is more difficult than high-precision relative mea- 
surements performed at the same time using the same 
instrument in the strong lensing case. 

Unfortunately, as Loeb (1998) pointed out, differences 
in the conditions of the images (e.g., shape, magnifica- 
tion, noise, etc.) of the same galaxy could already cause 
larger differences in the redshift measurements; even 
more problematically, proper motion of the lens could in- 
duce a difference in image redshifts (Birkinshaw & Gull 
1983) far exceeding that associated with the time de- 
lay. Therefore, even though the strong lensing case has a 
few advantages, it is still impractical to use the redshift 
difference and time delay between images of a strongly 
lensed source to measure the Hubble expansion at the 
source rcdshift. 
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APPENDIX 



MODEL STRUCTURES 

We require the curvature function K of the model structures to have a continuous second derivative with respect to 



r. For compensated models, we specify K directly 

K 



< 



Ai 






\ 3 r 


( ru 


-r \ 




-n J 



-2ri+7-n)(2r-7-i-ni) 
{ru-ri)' 2 



< r < n 
r\ < r < ru 
r > m 



(Al) 



where K\ = Si^af HfQ m , and (5i i? : is the initial overdensity of region I. For uncompensated models, we let the initial 
overdensity Si follow the same form as K (with K~i replaced by 8i t i) and then use Equations (4), (8), and (2) to obtain 
K . We choose ni for uncompensated models so that the uncompensated overdensity profiles roughly match those of 
compensated models in the central region. 

The left panel of Figure 4 shows overdensity profiles of the compensated (dotted lines) and uncompensated (solid 
lines) cluster models. For both models, the initial overdensity of region I at a, is set to 5n = 7.0 x 10 -4 , so that 
Si — 0.36 at redshift z = 0.5 to match the mean overdensity of overdense regions of the clusters found in Granett et al. 
(2008). The density of the uncompensated cluster grows with time everywhere. The void models behave similarly 
(not shown) but with inverted profiles and smaller amplitudes. Moreover, the uncompensated void develops slightly 
positive overdensity (S <C \S\\) at r > ru at late time. 

The right panel of Figure 4 illustrates the evolution of the function R/ar for the compensated cluster (solid lines) 
and the compensated void (dotted lines). The void has = —5.5 x 10~ 4 and S\ = —0.19 at z = 0.5, matching the 
mean overdensity of underdense regions of the voids in Granett et al. (2008). In our setup, underdense regions always 
expand faster than the background after t.j. Hence, for the voids, R(t,r) > ar at t > and the equality takes place 
where the mean overdensity within r vanishes. Similarly, for the clusters, R(t,r) < ar at t > U. 
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